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1. Introduction 

This is a paper of physical mathematics, which means that it addresses a problem of mathe- 
matics using tools of (theoretical) physics. The problem of mathematics is a venerable one, 
more than a hundred years old, namely the counting of (topologically inequivalent) knots. 
The physical tools are combinatorial methods developed in the framework of field theory 
and so-called matrix models. For a review of the history and recent developments of the 
first subject, see |jl|]. For an introduction for non physicists to matrix integral techniques, 
see for example [0-0. 

In this note, we show that by combining results obtained recently in knot theory and 
older ones on matrix integrals, and by using graphical decompositions familiar in field 
theory, one may reproduce and somewhat simplify the counting of alternating tangles and 
links performed in [Q. In section 2, we recall basic facts and definitions on knots and their 
planar projections; we also recall why integrals over large matrices are relevant for the 
counting of planar objects. 

Specifically, we shall consider the following integral 



over N xN complex matrices, in the large limit. In that limit, the integral is represented 
in terms of planar Feynman diagrams, with directed edges and four-valent vertices of the 
type which exhibit a close similarity with alternating knot diagrams in planar projec- 
tion, with crossings represented as /\. Thus the counting of planar Feynman diagrams 
(with adequate conditions and weights) must be related to the counting of alternating 
knots. A substantial part of this paper (section 3) is devoted to eliminating irrelevant or 
redundant contributions of Feynman diagrams. Once this is achieved, the results of 
are recovered. In the concluding section, we comment on the possible extensions of these 
methods. 

The observation that planar Feynman diagrams generated by matrix models can be 
associated to knot diagrams was already made in the matrix integral proposed there 
was more complicated, so that no explicit calculation was carried out. 




(1.1) 
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2. Basics 



2.1. Knots, links and tangles 

In this section, we briefly recaU some basic concepts of knot theory, referring to the hter- 
ature for more precise definitions. A knot is a smooth circle embedded in M^. A link is a 
coUection of intertwined knots (in the following, we shall not consider "unlinks" , i.e. links 
which can be divided in several non-intertwined pieces). Both kinds of objects are consid- 
ered up to homeomorphisms of JR^. Roughly speaking, a tangle is a knotted structure from 
which four strings emerge: if this object is contained in a ball B with the four endpoints of 
the strings attached on dB, topological equivalence is up to orientation preserving home- 
omorphisms of B that reduce to the identity on dB. The fundamental problem of knot 
theory is the classification of topologically inequivalent knots, links and tangles. 




Fig. 1: (a): a non prime link; (b): an irrelevant (or "nugatory") crossing; 
(c): 2 particle- reducible tangles, horizontal or vertical sums of two tangles 

It is common to represent such objects by their projection in a plane, with under/over- 
crossings at each double point and with the minimal number of such crossings. To avoid 
redundancies, we can concentrate on prime links and tangles, whose diagrams cannot be 
decomposed as a sum of components (Fig. 1) and on reduced diagrams that contain no 
irrelevant crossing. 

A diagram is called alternating if one meets alternatively under- and over-crossings 
as one travels along each loop. Starting with eight (resp six) crossings, there are knot 
(resp link) diagrams that cannot be drawn in an alternating way. Although alternate links 
(and tangles) constitute only a subclass (asymptotically subdominant), they are easier to 
characterize and thus to enumerate or to count. 




Fig. 2: The flype of a tangle 
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A major result conjectured by Tait and proved by Menasco and Thistlethwaite is that 
two alternating reduced knot or link diagrams on the sphere represent the same object if 
and only if they are related by a sequence of moves acting on tangles called "flypes" (see 
Fig. 2). 

The subproblem that we shall address here is thus the counting of alternating prime 
links and tangles. 

2.2. Matrix Integrals in the large N limit 

As a prototype of matrix integrals, we consider the integral over N x N hermitian matrices 

dMexp-A^tr Qm^ - |m4^ . (2.1) 

In order that the integral makes sense, the sign of g should originally be chosen negative. 
As is well-known, the large limit allows the analytic continuation to positive values 
of g. The series expansion of Z in powers of g may be represented diagrammatically by 
Feynman diagrams, made of undirected edges or "propagators" (2-point functions of the 
Gaussian model) with double lines expressing the conservation of indices, {MijMki)o = 

P "m 

-'k = jfSudjk, and of four-valent vertices ""y^ = gN5qi5jk5irnSnp- In the large N 



limit a counting of powers of N shows that the leading contribution to log Z is given by a 
sum of diagrams that may be drawn on the sphere , called "planar" by abuse of language. 
More precisely 

lim -^\ogZ= V weight g"" (2.2) 

planar diagrams 
with n vertices 

with a weight equal to one over the order of the automorphism group of the diagram. Once 

this has been realised, it is simpler to return to a notation with simple lines and rigid 

vertices It is this property of matrix integrals to generate (weighted) sums over planar 
diagrams that we shall use in the context of knot theory. 



3. From planar Feynman diagrams to links and tangles 

3.1. The matrix integral 

As mentionned in the Introduction, in the context of knot theory, it seems natural to 
consider the integral ( |1 . 1|) over complex (non hermitian) matrices, in order to distinguish 
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between under-crossings and over-crossings. However, this integral is closely related to the 
simpler integral (|2.1j ) in the large limit. Let us define the partition functions 



Z^^\a,g) = J dMdMUxp~Ntr(^aMM^ - ^{MM^f^ (3.1a) 
Z{a,g)= f dMexp-Ntr (^aM^ - ^mA (3.16) 



and the corresponding "free energies" 

F^fa a) - hm ^ log Z^'\a,g) 

The constant a can be absorbed in a rescaling M a~^M: 

Z(a,g) = a- — Z(l,-^) (3.3a) 
F(a,^7) = F(l,4) (3.36) 

and similarly for Z*-^-* and F^^\ However, it will be useful for our purposes to keep the 
parameter a. 

The Feynman rules of (|3.1a|) and those of ( p.l6|) are quite similar: the former have 



already been described in Sect. 2.2, while the latter are given by a propagator (MijM^^)o 



j ■=)>=. i = -j^SiiSjk and a four-vertex equal to gN times the usual delta functions. 

The only difference lies in the orientation of the propagator of the complex theory. However, 
in the large "planar" limit, this only results in a factor of 2 in the corresponding free 
energies, which accounts for the two possible overall orientations that may be given to the 
lines of each graph of the hermitian theory to transform it into a graph of the non-hermitian 
one. Therefore 

F^^Ha,g)=2F{a,g) (3.4) 

In addition to the partition functions and free energies, one is also interested in the 
correlation functions G'2n(a,^) = (^trM^^) and G^^^ia^g) = (^tr(MMt)'^). The first 
ones, namely the 2-point and 4-point functions, are simply expressed in terms of F. 

G4 = 4^ Gi'^=2^— (3.5) 

og ag 



so that G^l^ = G4 in the large N hmit. Furthermore, combining (^]^) with the homogeneity 



property ( p.3|) , one finds 

G2 = -- 2^F{a,g) = -(1 + gG^) (3.6a) 
a oa a 

G^i^ = -{l + gG^l^) (3.66) 
a 

which in particular proves that in the large limit G2 = G^2^ . 

In that "planar" limit, the function F has been computed by a variety of techniques: 
saddle point approximation [|^, orthogonal polynomials 0, "loop equations" (see for a 
review). With the current conventions and normalizations, 

F{a,g) = \ loga^ - ^(a^ - 1)(9 - a^) (3.7) 



where c? is the solution of 



3470^-0^ + 1 = (3.8) 



which is equal to 1 for (7 = 0. We have the expansion 



F{Xg) = \[g^\g'^%^+^^g'^...) 

00 



- y (3^)^^^^ (3.9) 



The 2- and 4-point functions are thus 



G2 = G^ = — a2(4 - a^) (3.10a) 
3a 



^4 = ^a4(3 - a^) (3.106) 



Gi = G4 - 2Gl = \.a^(a^ - l)(2a^ - 5) (3.10c) 

where G^ is the connected 4-point function. More generally, G2n is of the form a~'^ 
times a polynomial in a^. The singularity of a'^ at g/a'^ = 1/12 determines the radius of 
convergence of F and of the G2n- 



3.2. Removal of self- energies 

If we want to match Feynman diagrams contributing to the free energy with prime 
knots or hnks, we have to ehminate a certain number of redundancies. We have first 
to ehminate the "self-energy insertions" that correspond either to non prime knots or 
to irrelevant crossings ("nugatory" in the language of Tait). This is simply achieved by 
choosing a as a function of g such that 



G2{oi{g),g) = 1 . 



(3.11) 



The function a^i^g) := a?{a{g),g) is obtained by eliminating a between Eqs. ( p.8| ) and 



S^ra^ - + 1 = 



1 



a^(4 - a^) = a 



(EH' 



which implies that o?{g) is the solution of 



21 g = (a^ - l)(4-a' 



2\2 



(3.12) 



equal to 1 when g = 0; «((?) is then given by 



aig) = ^a\g)i4-a\g)) 



(3.13) 



(a) 



+ — 



a a a a a 



(b) 





9 G 



9G 



Fig. 3: (a): decomposition of the two-point function into its one-particle- 
irreducible part; (b) discarding one- vertex-reducible contributions 
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Let us now consider correlation functions. In field theory, it is common practice to 
define truncated diagrams, whose external lines carry no propagator, and one-partide- 
irreducible ones, that remain connected upon cutting of any line. 

The two-point function G2{o(,g) may be expressed in terms of the "self-energy" func- 
tion E 

G2{a,9) = r (3.14) 

a — zj(q;, g) 

which is the sum of (non-trivial) truncated, one-particle-irreducible graphs (Fig. 3). One 
can further discard the contributions that are one-vertex-reducible by defining E' (see 
Fig. 3(b)) 

E'{a,g) = E{a,g)-2G2{a,g) (3.15) 

If we now remove the self-energy insertions by imposing condition ( |3.11| ) , we find that 
Eqs. ( p.l4| ) and ( |3.15| ) simplify, so that T,{g) := 'E{a{g),g) and ^'(^f) :— T^' {a{g), g) are 
given by 

E{g) = a{g) - 1 (3.16a) 
E'{g)=a{g)-l-2g (3.166) 

i.e. obtained from a{g) by removing the first terms in its expansion in powers of g. 

The procedure extends to all correlation functions. Finally one obtains the corre- 
sponding "free energy" F^^\g) by dividing the term of order n of T,'{g) by 2n (= number 
of times of picking a propagator in a diagram of the free energy to open it to a two-point 
function, cf Eq. ( |3.6|) above) . 

We may also compute the function r{a,g) = 02(0:: g)~^G^{a{g), g), which counts 
the truncated (automatically one-particle-irreducible) connected 4-point functions. After 
removal of the self-energy insertions, T{g) := r{a{g),g) becomes simply (Eqs. (|3.6a|) , 
( P^) and ( p^) ) 



Explicitly, 



Tig) = ^=2^F(^\g) (3.17) 
9 dg 



m = - (4_j2(^))2 («'(^) - l)(2a'(^7) - 5) (3.18) 
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Perturbatively one finds 



a 



;i, g) = l + 3g + 18g'^ + 135g^ + 113^ + 10206£?^ + 96228/ + ■ 



^2(1,^7) 



(a) 



(b) 



i 




l + 2g + 9g^ + 54g^ + 378g^ + 2916£?^ + 24057£?^ + 
1 + 2^ + / + 2^3 + 6/ + 22^^ + 91/ + ■ • • 
g + 2g^ + 6/ + 22/ + 91^^^ + 408/ + ■ ■ • 



4 3 4 



+ 



11/ , 91g 



+ 



12 



+ ■ 








5^ 



(3.19) 



(c) 



X >cx 



x:>c:x 



Fig. 4: (a): the first links, with the labeUing of |jT^; (b) the corresponding 
diagrams contributing to F^^\ For simphcity, the diagrams are not oriented, 
but the weights are those of the (MM^)"^ theory; (c) diagrams up to order 3 
contributing to F: the last four are pairwise flype-equivalent. 



The first terms of F^^^ match the counting of the first prime links weighted by their 
symmetry factor (see Fig. 4). 

F^'Hg) = -/ + + {\ + \)9' + + 1 + ^)9' + (3.20) 

Starting with order 6, however, there is an overcounting of links due to neglecting the flype 
equivalence. The overcounting occurs already at order 3 if F is used to count tangles. 

Asymptotic behavior of the coefficients fn of F^^\g) 
The singularity is now given by the closest zero of the equation g/a^[g) = 1/12, which 
gives g^ = ^. Thus 

~ const 6"n-5 (3.21) 

with 

97 

6.75 . (3.22) 
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3.3. Two-particle irreducibility 

Since the flype acts on tangles, we have to examine more closely the generating function 
r{g) of connected 4-point functions with no self-energies. We want to regard the cor- 
responding diagrams as resulting from the dressing of more fundamental objects. This 
follows a pattern familiar in field theory, whose language we shall follow, while indicating 
in brackets the corresponding terminology of knot theory. 

We say that a 4-leg diagram (a tangle) is two-particle-irreducible (2PI) (resp two- 
particle reducible, 2PR) if cutting any two distinct propagators leaves it connected (resp 
makes it disconnected). A 2PR diagram is thus the "sum" of smaller components (see 
Fig. 1). A fully two-particle-irreducible diagram is a 4-leg diagram such that any of its 4- 
leg subdiagrams (including itself) is 2PI. Conversely a fully 2PR diagram (algebraic tangle) 
is constructed by iterated sums starting from the simple vertex. We shall also make use 
of skeletons, which arc generalized diagrams (or "templates" ) in which some or all vertices 
are replaced by blobs (or "slots"). The concepts of fully 2PI skeleton (or basic polyedral 
template) and of fully 2PR skeleton (algebraic template) follow naturally, with however 
the extra condition that only blobs should appear in the former. 

These blobs may then be substituted by 4-leg diagrams, resulting in a "dressing". As 
will appear, the general 4-leg diagram (tangle) results either from the dressing of a fully 
2PI skeleton by generic 4-leg diagrams (type I tangles) or from the dressing of non trivial 
fully 2PR skeletons by type I tangles. The action of fiypes will be only on the fully 2PR 
skeletons appearing in this iteration. 

Because in the {MM'^)^ theory with 4-valent vertices, there is no diagram which is 
two-particle-reducible in both channels, any diagram of T'{g) must be 

i) 2-particle-irreducible in the vertical channel (V-2PI) but possibly 2-particle-reducible 
in the horizontal one. We denote V{g) the generating function of those diagrams. 

ii) or 2-particle-irreducible in the horizontal channel (H-2PI) but possibly 2-particle- 
reducible in the vertical one. We denote H{g) the generating function of those dia- 
grams. 

Obviously there is some overlap between these two classes, corresponding to diagrams that 
are 2-particle-irreducible in both channels (2PI), including the simple vertex. Let D{g) 
denote their generating function. We thus have 

T = H-\-V-D. (3.23) 
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Now since V encompasses diagrams that are 2PI in both channels, plus diagrams that 
are once 2PR in the horizontal channel (i.e. made of two H-2PI blobs joined by two 
propagators), plus diagrams twice H-2PR, etc, we have (see Fig. 5) 



V = D + HH + HHH + ■■■ 

= D+ . (3.24) 










+ 






+ ( 




\ 











Fig. 5: Graphical representation of Eqs. ( |3.23| ) and ( p.24| ). 
Since for obvious symmetry reasons, H = V , v^e have the pair of equations 

T = 2H -D 



H = D + 



Eliminating H yields 



1-H 
1 -F 



1 + F 



(3.25) 



(3.26) 



Eq. (|3.26|) can be inverted to allow to reconstruct the function r{g) out of the 2PI function 
D{g). For later purposes, we want to distinguish between the single vertex diagram and 
non-trivial diagrams: 

D(g)=g + C{9) (3.27) 



In terms of these new variables, one has: 

ng)=r{gXi9)} 

(note the braces are used to distinguish the different variables) where 

1 



T{g, C} = - [l - ^7 - C - V(l-^7-C)'-% + C) 



is the generating function 



(3.28) 



(3.29) 



(3.30) 
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of the number of fully 2PR skeleton diagrams with m vertices and n blobs (algebraic 
templates, including the trivial one made of a single blob and no vertex). Similarly, one 
could define H{gX} = VigX}^ ^ind of course D{gX} = 9 + C- Note that the function 
T{g, C} does not depend on the precise form of r{g), since we have only used the relation 
( pi.26| ) which was derived from general diagrammatic considerations. 

Inversely, we shall need the fully 2PI skeletons which are obtained from ({g) by re- 
placing subgraphs that have four legs with blobs. Defining the inverse function g[r] of 
T{g), the generating function of these skeleton diagrams is simply C[r] := C((7[r]), or more 
explicitly 



C[r] 



,1 



i + r 



The function g[r] satisfies by definition r(5f[r]) 



between Eqs. (11) and (16). Setting 77 = 1 — a , we recover the system of 



(3.31) 

r, and is found by eliminating 



which leads to 



and finally 



9[r] = o 



27g 

r 

1 1 



-77(3 + 77) 
1 



(3.32) 



77(3 + 277) 



2 (r + 2)3 



C[r] 



1 + lor - 2r^ - (1 - 4r) 



+ 2-r-^7(r) 



i + r 



(3.33) 

(3.34) 

this is 



is the desired generating function of fully 2PI skeletons (in the notations of | 
q{g))- The property mentionned above that F is obtained by dressing is expressed by the 
identity 

r{g) = am] + ^rn,n9'^c[m] = r{^7, a9)} . (3.35) 



Perturbatively, we find 

T{g) = g + 2g^ + 6g^ + 22/ + 91^7^ + ■ ■ ■ 

D{g) = g + g' + lOg"" + 74^7^ + 492(7^ + ■ ■ ■ 

r{^7,C} = ^7 + C + 2(^7 + 0' +6(^7 + 0' + 22(^7 + 0^90(^7 + C)' + --- 
g[T] = T-2T^ + 2T^ - 2T^ + - 2F^ - 2F^ - 8F^ - 22F^ - 68F^° 



(3.36) 



+ 



C[F] = F^ + 4F^ + 6F^ + 24F9 + 66F^° + ■ ■ ■ . 
The first contributions to ([T] are depicted on Fig. 6. The closest singularity of g[r] and 
of C[F] is F* = 1/4 which corresponds to ^^[F*] = g^ = 4/27 and CP*] = 1/540. 



11 



Fig. 6: The first contributions to C[r]. All diagrams can be obtained from 
tlie ones depicted by rotations of 90 degrees. 



3.4- Quotienting by the flype 

The last step is to take into account the flype equivalence; it is borrowed from the discussion 
of Sundberg and Thistlethwaite and reproduced here for completeness. The fully 2PR 
skeletons contained in Eq. ( |3.29D now have to be replaced with skeletons in which the 
flype equivalence has been quotiented. Then, they have to be dressed by 4-point 2PI 
functions using Eq. ( p.34|) . 

Let r{(7, C} be the generating function of the number ^rnn of flype- equivalence classes 
of fully 2PR skeleton diagrams (algebraic templates) with m vertices and n blobs 

^{9X} = Y.'im,r.9'^C ■ (3.37) 

Let H{g, C} (resp. V{g, C}) denote the generating function of flype-equivalence classes of 
skeletons which are 2PI in the horizontal (resp. vertical) channel, including the single blob 
and the single vertex. In a way similar to the decomposition performed above for F, cf 
Eq. (|3.23| ), we write 

t = H + V-D (3.38) 



where D{g, C} = fi* + C- The equation analogous to Eq. (|3.24| ) is 



V = D + gt + {H-gf + {H-gf + --- 

^D + gr + — (3.39) 

where by flyping we can remove the simple vertices inside the H and put them as a single 
contribution gV. 

As before, H = V for symmetry reasons, and after eliminating it one gets an algebraic 
equation for r{g, C} 

f2 - (1 + ^7 - Of + C + = (3.40) 
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with solution 



[l + g-C)-J(l-g + Cr-8C^ ° 



1-9 



(3.41) 



which should be compared with Eq. ( p.29|) . 

The last step to get the generating function 

tig):=t{gra9)} (3.42) 

of flype equivalence classes of (alternating, prime) tangles is to define ({g) via the relation 
(PI), i.e. 

Cig) = am] (3.43) 

(The reason we can use Eq. ( p.34|) without any modification is that fully 2PI skeleton 
diagrams are not affected by flypes.) In other words, T{g) is the solution to the implicit 
equation 

t{g)^f{gX[r{9)]} , (3.44) 
where ([T] and r{gX} are provided by Eqs. (|3.33| ), ( |3.34| ) and ( p.41|) , and which vanishes 
aX g = 0. 

Perturbatively 

f(^) =g + 2g^ + Ag^ + IQg^ + 29^7^ + 98^7^ + 372(7^ + ■ • ■ . (3.45) 

Asymptotic behavior of T[g): we know the closest singularity of C[r]; we therefore set 
f = = 1/4 and C = C* = 1/540 in Eq. (p^) and solve for g: 



^ -101 + 721001 

270 ^ ' 

This provides the asymptotic behavior of the coefficients /„ of the "free energy" F^^\g) 
defined by f(^7)=2^F(i)(^7): 

fn ~ const n~i (3.47a) 



h = (101 + V21001)/40 ^ 6.147930 . (3.476) 

Note that while the first factor 6" in the expansion ( "bulk free energy" ) is non-universal, 
the second factor n~^/^ is universal (critical behavior of pure gravity [Q; the exponent 
ct = — I gives the "string susceptibility" exponent 7 of non-critical string theory ||12|| : 
7 = a-|-3 = — |) and in particular is identical in Eqs. (|3.21|) and (|3.47a|) . Note also that 



this free energy counts links with a weight that takes into account the equivalence under 
flypes and which is less easy to characterize. 
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4. Concluding remarks 

In this paper we have reproduced the results of 0] using prior knowledge of graph counting 
derived from matrix models, while Sundberg and Thistlethwaite were using results of Tutte 



TlJ]. Admittedly the progress is modest. We hope, however, that our method may give 



some clues on problems that are still open, such as the counting of knots rather than 
links. To control the number of connected components of a link is in principle easy in our 
approach. We should consider an integral over n matrices Mq, , a = 1 , ■ ■ ■ , n interacting 
through a term ^(M^M^)^ and look at the dependence of F^^^^ on n. The term linear in 
n receives contributions only from one-component diagrams, hence after a treatment similar 
to that of Sect. 3, it should give a generating function of the number of knots, weighted as 
before by their symmetry factor. Unfortunately, the computation of these matrix integrals 
and their subsequent treatment (removal of self-energies and flype equivalences) is for 
generic n beyond our capabilities (see however []T^ for a first step in this direction). 

Another problem on which matrix technology might prove useful would be in the 
counting of non alternating diagrams. But there the main problem is knot-theoretic rather 
than combinatorial: how does one get rid of multiple counting associated with Reidemeister 
moves? 
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